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Abstract. In recent years, much studies have focused on 
determining the origin of the large-scale line-profile and/or 
photometric patterns of variability displayed by some ap- 
parently single Wolf-Rayet stars, with the existence of an 
unseen (collapsed?) companion or of spatially extended 
wind structures as potential candidates. We present ob- 
servations of WR 1 which highlight the unusual character 
of the variations in this object. Our narrowband photo- 
metric observations reveal a gradual increase of the stellar 
continuum flux amounting to Av ~ 0.09 mag followed by 
a decline on about the same timescale (3-4 days). Only 
marginal evidence for variability is found during the 11 
following nights. Strong, daily line-profile variations are 
also observed but they cannot be easily linked to the pho- 
tometric variations. Similarly to the continuum flux vari- 
ations, coherent time-dependent changes are observed in 
1996 in the centroid, equivalent width, and skewness of 
He II A4686. Despite the generally coherent nature of the 
variations, we do not find evidence in our data for the 
periods claimed in previous studies. While the issue of a 
cyclical pattern of variability in WR 1 is still controversial, 
it is clear that this object might constitute in the future 
a cornerstone for our understanding of the mechanisms 
leading to the formation of largely anisotropic outfiows in 
Wolf-Rayet stars. 
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1. Introduction 

It is now fairly well-established that apparently single 
Wolf-Rayet (WR) stars may display two distinct (but 
probably non mutually exclusive) spectroscopic patterns 
of variability: (a) small-scale emission features systemat- 
ically moving from the line center to the line wings on 
an hourly timescale (e.g., Lepine 1998[ ); (b) dramatically 
larger line-profile deformations operating on a much longer 
basis (~ days, e.g., [Smith fc Willis 1994| ). Ahhough the 
first phenomenon, observed in most (if not all) WR stars, 
is believed to be the consequence of the fragmented, possi- 
bly turbulent nature of the outflow, the origin of the latter 
type of variability is still very much elusive. 

Remarkable in this respect, is the existence of a well- 
established (although strongly epoch-dependent) large- 
scale, cyclical pattern of variability in the two appar- 
ently single WR s tars WR 6 {V = 3.763 ± 0.002 d; 
Firmani et al. 1980D and WR 134 {V = 2 .27 ± 0.04 
d; [McCandliss et al. 1994^ [Morel et al. 19991). A major 



observational effort has been directed on establishing 
the true nature of these peculiar objects, i.e., whether 
this cyclical variability is induced by an orbiting un- 
seen (collapsed?) companion or by the rotational modula- 



tion of large-scale wind structures (e.g., Vreux et al. 1992 



Morel 1998; Harries et al. 199!;, and references therein) 



Although the exact nature of these stars has yet to be 
unambiguously settled, these studies reveal that rota- 
tional modulation constitutes an attractive alternative to 
the binary hypothesis, especially considering the recent 
recognition that some O stars (the progenitors of WR 
s tars) might possess s uch azimuthally structured outflows 
( [FuUerton et al. 1997[ ; [Kaper et al. 1997[ ). 

A prime target for further investigations is the seldom- 
studied, apparently single WN 5 star WR 1 (HD 4004) 
that has recently been shown to present a spectral 
pattern of variability very similar to that of WR 6 
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( [Niedzielski 1995|, 199 6a, b; Wessolowski & Niedzielski 
1996; Niedzielski 1998). Strong line-profile variability was 
observed, as well as apparently cyclical (according to V 
« 2.667 days) variations in the EWs of Heii A4686 and 
He II A5412 ( [Niedzielski 1996aD . The first claim of peri- 
odic variability in WR 1 with V ~ 7.7 days was made by 



Lamontagne (1983) from an analysis of the radial veloc- 
ity variations of Heii A4686. The first photometric mon- 
itoring of this object has shown WR 1 to be variable, 
with an indicati on of a 6.1 day period (M offat & Shara 
1986). Recently, Marchenko et al. (1998a) discussed Hip- 
parcos broadband photometric data which revealed that 
WR 1 also displays relatively long-term photometric vari- 
ations, with a marginal evidence for a. V ^ 11.68 ± 0.14 
day periodicity. 

As can be seen, controversy persits in the literature 
concerning the possible cyclical nature of the variations in 
WR 1. We present in this paper the results of spectroscopic 
and photometric monitoring of WR 1 carried out in 1995 
and 1996 aiming at shedding some light on this issue. 

2. Observations and Reduction Procedure 

8.1. Photometry 

The photometric variability of WR 1 has been inves- 
tigated during the interval 1996 September 18-October 
5 by use of the single channel photometer Cuentapul- 
sos on the 0.84 m telescope of the Observatorio As- 
tronomico National at San Pedro Martir (Mexico). Two 
additional objects were monitored during this observing 
run, namely, WR 3 and WR 153. The nights were gen- 
erally clear. WR 1 was observed through a narrowband 
V filter centered on 5140 A (FWHM = 90 A). This filter 
samples a continuum-dominated region of the WR spec- 
trum. We applied the following sequence of 60 s integra- 
tion through a 25" diaphragm: sky, C2, CI, WR, CI, WR, 
CI, C2, sky. The same nearby comparison stars as used 



by Moffat & Shara (1986) have been chosen. These com- 
parison stars are similar in terms of their magnitude and 
colour to WR 1: AB (WR - CI) = - 0.17, A[B ~ V] (WR 

- CI) = - 0.17, AB (WR - C2) = - 0.26, A[B ~ V] (WR 

- C2) = - 0.19. An extinction coefficient — 0.20 was 
used throughout the data reduction. The scatter in the 
(C2 — CI) data for the whole dataset amounts to cr = 4.7 
mmag. The differential magnitudes quoted in Table 1 are 
averaged over two consecutive cycles typically separated 
by about 20 minutes. 

2.2. Spectroscopy 

2.2.1. Long-slit and Reticon Spectra 

Long-slit spectra of WR 1 have been obtained during vari- 
ous campaigns at the Observatoire du Mont Megantic and 
Dominion Astrophysical Observatory (Canada) in 1995 
October and 1996 September. Reticon spectra were also 



obtained at DAO in 1996 November. The 1996 campaign 
at the Observatoire du Mont Megantic was coordinated 
to support the photometric campaign described above. 
Table 2 lists the mode of observation, the dates of the 
spectroscopic observations, the interval of the observa- 
tions in heliocentric Julian dates, the observatory name, 
the number of CCD spectra obtained, the selected spec- 
tral domain, the reciprocal dispersion of the spectra, and 
the typical signal-to- noise ratio (S/N) in the continuum. 
The spectra were reduced using the IRAf]^ data reduction 
packages. The bias and sky subtraction, flat-field division, 
removal of cosmic ray events, extraction of the spectra, 
and wavelength calibration were carried out in the usual 
way. Spectra of calibration lamps were taken immediately 
before and after the stellar exposure. The stellar spec- 
tra were subsequently continuum normalized by fitting a 
low-order Legendre polynomial to carefully selected line- 
free regions. In order to minimize the spurious velocity 
shifts induced by an inevitably imperfect wavelength cal- 
ibration, the spectra were coaligned in velocity space by 
using the interstellar doublet Nai AA5890, 5896 as fidu- 
cial marks. When not available, the doublet Can AA3934, 
3968 or the diffuse interstellar band at 4501 A were used. 
A trend for a systematic shift of the zero point of the 
wavelength scale has been corrected by redshifting most 
of the spectra by an average value of 35 km s~^. 

2.2.2. Echelle Spectra 

Echelle spectra have been obtained during the period 1996 
September 16-19 with the Echelle spectrograph Espresso 



(Levine & Chakrabarty 199£) on the 2.1 m telescope of 
the San Pedro Martir Observatory. The UCL camera and 
a 1024 X 1024 coated CCD-Tek chip have been used. The 
selected grating (300 lines mm"-'^) yields a reciprocal dis- 
persion of 0.16 and 0.23 A pixel~^ at H7 and Ha, respec- 
tively. The spectra cover 27 orders and span the spectral 
range 3720-6900 A. 

The reduction procedure (bias subtraction, division by 
a normalized flat field, removing of scattered light, ex- 
traction of the orders) was carried out using the reduction 
tasks in the IRAF package echelle. Comparison spectra of 
Th-Ar lamps have been used for the wavelength calibra- 
tion. The typical accuracy of the wavelength calibration 
can be judged by the dispersion in the heliocentric radial 
velocities of the interstellar line Nai A5890: cr w 2 km s^^. 
The instrumental response has been removed by fitting a 
low-order Legendre polynomial to the continuum sections 
of each order of an 0-star spectrum, and then dividing the 
corresponding order of the WR 1 spectra by this polyno- 
mial. The resultant, almost flat orders were then combined 



^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, operated by the Association of Universities for Re- 
search in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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HJD 
(- 2,449,000) 


WR 


- CI 


WR 


- C2 


C2 - 


- CI 


HJD 

(- 2,449,000) 


WR 


- CI 


WR 


- C2 


C2 - CI 


1345.721 


- 


043 


- 


114 


+ 





071 


1356.693 


- 


054 


- 


121 


+ 


0.067 


1347.702 


- 


074 


- 


134 


+ 





061 


1356.773 


- 


057 


- 


128 


+ 


0.071 


1347.731 


- 


063 


- 


121 


+ 





058 


1356.847 


- 


065 


- 


130 


+ 


0.065 


1347.760 


- 


064 


- 


137 


+ 





073 


1356.902 


- 


054 


- 


127 


+ 


0.073 


1347.846 


- 


065 


- 


134 


+ 





069 


1356.984 


- 


054 


- 


130 


+ 


0.076 


1347.882 


- 


066 


- 


141 


+ 





075 


1357.656 


- 


064 


- 


139 


+ 


0.075 


1348.807 


- 


095 


- 


163 


+ 





068 


1357.723 


- 


064 


- 


129 


+ 


0.065 


1348.831 


- 


105 


- 


173 


+ 





067 


1357.781 


- 


055 


- 


121 


+ 


0.066 


1348.857 


- 


094 


- 


170 


+ 





076 


1357.847 


- 


058 


- 


127 


+ 


0.070 


1348.886 


- 


103 


- 


161 


+ 





058 


1357.914 


- 


057 


- 


129 


+ 


0.072 


1349.678 


- 


126 


- 


197 


+ 





071 


1358.735 


- 


055 


- 


126 


+ 


0.071 


1349.761 


- 


119 


- 


185 


+ 





066 


1358.773 


- 


059 


- 


125 


+ 


0.066 


1349.854 


- 


129 


- 


195 


+ 





065 


1358.829 


- 


048 


- 


109 


+ 


0.061 


1349.933 


- 


127 


- 


193 


+ 





065 


1358.881 


- 


061 


- 


123 


+ 


0.062 


1349.991 


- 


130 


- 


194 


+ 





064 


1358.943 


- 


051 


- 


117 


+ 


0.066 


1350.873 


- 


101 


- 


164 


+ 





063 


1358.992 


- 


058 


- 


120 


+ 


0.062 


1350.947 


- 


104 


- 


167 


+ 





063 


1359.642 


- 


051 


- 


118 


+ 


0.067 


1351.006 


- 


100 


- 


164 


+ 





064 


1359.702 


- 


040 


- 


112 


+ 


0.072 


1352.851 


- 


059 


- 


123 


+ 





064 


1359.767 


- 


043 


- 


115 


+ 


0.071 


1352.913 


- 


053 


- 


126 


+ 





073 


1359.832 


- 


049 


- 


116 


+ 


0.067 


1352.953 


- 


054 


- 


119 


+ 





065 


1359.909 


- 


051 


- 


117 


+ 


0.066 


1352.964 


- 


052 


- 


118 


+ 





066 


1359.958 


- 


050 


- 


122 


+ 


0.072 


1353.000 


- 


058 


- 


125 


+ 





067 


1360.004 


- 


059 


- 


122 


+ 


0.063 


1353.649 


- 


062 


- 


125 


+ 





063 


1360.632 


- 


063 


- 


133 


+ 


0.070 


1353.747 


- 


058 


- 


126 


+ 





067 


1360.697 


- 


046 


- 


114 


+ 


0.067 


1353.836 


- 


048 


- 


115 


+ 





067 


1360.761 


- 


059 


- 


126 


+ 


0.067 


1353.936 


- 


049 


- 


115 


+ 





066 


1360.842 


- 


057 


- 


116 


+ 


0.060 


1354.013 


- 


060 


- 


126 


+ 





066 


1360.912 


- 


061 


- 


123 


+ 


0.062 


1354.671 


- 


051 


- 


124 


+ 





073 


1361.637 


- 


060 


- 


122 


+ 


0.061 


1354.754 


- 


057 


- 


133 


+ 





076 


1361.694 


- 


060 


- 


132 


+ 


0.073 


1354.810 


- 


048 


- 


124 


+ 





076 


1361.768 


- 


061 


- 


121 


+ 


0.060 


1354.875 


- 


047 


- 


113 


+ 





066 


1361.851 


- 


065 


- 


124 


+ 


0.060 


1354.948 


- 


062 


- 


129 


+ 





067 


1361.952 


- 


065 


- 


131 


+ 


0.066 


1355.002 


- 


054 


- 


121 


+ 





068 


1361.997 


- 


058 


- 


126 


+ 


0.069 


1355.668 


- 


044 


- 


113 


+ 





070 


1362.625 


- 


053 


- 


118 


+ 


0.066 


1355.735 


- 


056 


- 


124 


+ 





068 


1362.720 


- 


060 


- 


124 


+ 


0.064 


1355.792 


- 


058 


- 


121 


+ 





063 


1362.846 


- 


052 


- 


121 


+ 


0.069 


1355.860 


- 


060 


- 


123 


+ 





063 


1362.910 


- 


055 


- 


123 


+ 


0.068 


1355.912 


- 


050 


- 


121 


+ 





071 


1362.983 


- 


053 


- 


123 


+ 


0.070 


1355.993 


- 


064 


- 


116 


+ 





052 

















Table 2. Journal of spectroscopic observations 



Mode of 




HJD 




Number 


Spectral 


Reciprocal Dispersion 




Observation 


Date 


(- 2,449,000) 


Observatory'^ 


of Spectra 


Coverage (A) 


(Apix-i) 


S/N 


Long-slit 


95 Oct 


992-1005 


OMM, DAO 


7 


3830-4300 


1.6 


« 160 










15 


4360-5080 


1.6 


« 170 










11 


5050-5945 


1.6 


« 290 


Echelle 


96 Sep 


1342-1345 


SPM 


72 


3720-6900 


0.16 (H-,) — 0.23 (H^) 


» 25 (H-,) — 50 (H„) 


Long-slit 


96 Sep 


1346-1352 


OMM 


13 


4400-5080 


1.3 


120 










14 


5050-6500 


1.3 


190 


Reticon 


96 Nov 


1386-1389 


DAO 


12 


5100-6110 


0.5 


« 75 



^ OMM: Observatoire du Mont Megantic 1.6 m; DAO: Dominion Astrophysical Observatory 1.8 m; SPM: San Pedro Martir Observatory 2.1 m. 



using the IRAF task scomb. This procedure proved to be 
satisfactory, except when the overlapping regions fall in 
a spectral domain with very steep intensity gradients, as 
in the blue wing of Hen A4686 (see FigJ|). The combined 
spectra were then subdivided into spectral regions roughly 
corresponding to those of the long-slit spectra discussed 
above (§2.2.1). For consistency purposes, these spectra 
have been continuum normalized by fitting a low-order 
Legendre polynomial to the same continuum sections that 
selected for the long-slit spectra. 



3. Results 



3.1. Photometric Variations 

The light curve of WR 1 is plotted as a function of the 
heliocentric Julian date of observation in Figure g. The 
main feature of this light curve is the gradual increase 
of the stellar continuum flux amounting to Aw « 0.09 
mag beginning at HJD 2,450,346, followed by its decline 
on about the same timescale after HJD 2,450,350. During 
the last 11 nights (after HJD 2,450,353), only marginal 
evidence for variability is found, with cr(WR — C)/a{G2 
— CI) « 1.3. This can be compared with a value of 



1.9 derived by Moffat fc Shara (1986) on the basis of 14 



nights of broadband B observations in 1984. Noticeable 
is the striking lack of recurrence in the continuum flux 
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Fig. 1. Two upper panels: differential v magnitudes of 
WR 1 in 1996 September relative to the two compari- 
son stars CI and C2. Lower panel: differential (C2 — CI) 
magnitudes. Dates when simultaneous spectroscopic ob- 
servations were performed are indicated by vertical tick 
lines in the uppermost panel. 
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Radial velocity (km s ^) 



Fig. 2. Superposition of the 21 EchcUc spectra acquired 
on 1996 September 18 (around HJD 2,450,344.873) for 
the spectral range encompassing Heii A4686. These spec- 
tra have been obtained on a total time interval of about 
5.5 hr. For comparison purposes, the dashed line shows 
the mean spectrum of the 1996 September dataset. The 
discontinuities observed at about — 2500 km and — 
1200 km s^^ are artefacts of an imperfect connection of 
the orders. The projected velocities, as everywhere in this 
paper, are heliocentric and refer to the line laboratory rest 
wavelength. 



data. A periodicity search in the (WR — CI) and (WR 
— C2) data acquired after HJD 2,450,353 yields no evi- 
dence for any periodic signals in the range 6-8 days that 
would possibly account for the variations observed be- 
fore this date. The highest peaks in the power spectra 
of these two datasets (after corre ction by the CLEAN al- 
gorithm; see Roberts et al. 1987 ) appear (with consider- 
able uncertainty) for periods of about 4 days. However, 
the significance of these periodic signals is very low and 
they must thus be regarded as spurious. The periods pro- 
posed by Niedzielski (1996a; 2.667 days) and Marchenko 
et al. (1998a; 11.68 days) are inconsistent with the global 
light-curve morphology. On the other hand, the periods 
claimed by Lamontagne (1983; 7.7 days) and Moffat & 
Shara (1986; 6.1 days) are only consistent with part of 
the data (before HJD 2,450,353). We will come back to 
this point later (§4). 

3. 2. An Overview of the Line-profile Variations 

In the following, we will generally illustrate the line-profile 
variations of WR 1 by means of the strong He ii A4686 line 



as its variations are qualitatively similar to those affecting 
the other Heii features (see §3.5; Niedzielski 1996a). One 
of the most outstanding features of the line-profile vari- 
ability is the relatively low level of variability observed 
within one night |^ Such a gradual pattern of variability 
has been observed whenever several spectra have been 
obtained during a single night of observation. This phe- 
nomenon is illustrated in Figure ^ where we show a super- 
position of the 21 Hell A4686 line profiles obtained on the 
night of 1996 September 18 (around HJD 2,450,344.873). 
In sharp contrast, however, noticeable night-to-night vari- 
ations are observed. This is illustrated in Figure y| which 
shows a montage of the nightly mean spectra obtained 
in 1995 October {left panel) and 1996 September {right 
panel) for the spectral region encompassing Hen A4686 
(it has to be kept in mind that because of the paucity of 



^ Note that although probably present in WR 1, the small- 
scale emission-excess features travelling on an hourly timescale 
on the top of the line profiles (probably induced in WR 
stars by outwardly moving shocked and/or turbulent material; 



Lepine & Moffat 1999) cannot be studied here, owing to the 
insufficient S/N achieved for the Echelle spectra. 
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-2000 2000 -2000 2000 

Radial velocity (km s^^) 

Fig. 3. Montage of the continuum-normalized, nightly mean spectra obtained in 1995 October {left panel) and 1996 
September {right panel) for the spectral range encompassing Hen A4686 (the number of individual spectra averaged 
to form the corresponding nightly mean is indicated between the two panels) . The mean heliocentric Julian date of the 
observations (- 2,449,000) is indicated on both sides of the panels. For comparison purposes, the unweighted means for 
1995 October (7 spectra) and 1996 September (8 spectra) are overplotted as dashed lines. The spectra of consecutive 
nights are offset by 2.5 units of continuum in the intensity scale. 



the photometric data, the strength of the spectral lines 
in these continuum-normalized spectra has not been cor- 
rected for the varying continuum flux level). 



3.3. Temporal Variance Spectrum Analysis 
The Temporal 



FuUerton ct al. 1996) 



Variance 
has been 



Spectrum (TVS; 
used to assess the 



level of spectral variability as a function of wavelength. 
The (square root of the) TVS, giving the typical "size" of 



the deviations from a template-weighted mean spectrum 
(expressed in percentage of the continuum flux) is shown, 
along with the 99.0 % confidence level for variability, in 
Figure ^. As can be seen, highly significant variability 
affects the emission lines, with the notable exception of 
N V A 4945 for which only little evidence for variability is 
found. Many peaks (i.e., locations of "preferential" vari- 
ability) can be distinguished in the TVS. The locations 
of these peaks (in terms of the projected velocity referred 
to the laboratory rest wavelength of the spectral feature 
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4600 4700 4800 4900 5400 5600 5800 

Wavelength (A) Wavelength (A) 

Fig. 4. Square root of the Temporal Variance Spectrum {short- dashed line), along with the threshold indicating the 
99.0 % confidence level for significant variability {long-dashed line). The mean of the spectra used in the calculations 
is overplotted in arbitrary units {solid line). 



Table 3. Projected velocities of the TVS subpeaks ob- 
served in He II A4686, Heii A4859, and Hen A5412.'* 



He 11 A4686 


Hen A4859 


Hen A5412 


- 2260 


- 1790 


- 1880 


- 1060 


- 1060 


- 1040 




- 550 




430 


-\- 290 


-)- 690 


+ 1560 


+ 1390 


+ 1430 


+ 2460 







" In km s ^ and referred to the line laboratory rest wavelength. 



in question) are quoted in Table 3. This investigation 
shows that: (a) the variability often extends to velocities 
comparable to the wind terminal velocity (uqo ~ 2135 
km s^^; Rochowicz fc Niedziclski 1995 ), e.g., t; w - 1790 
km s~^ for the relatively unblended line Heii A4859; (b) 
the TVS structure — made up of several subpeaks — 
presents some similarities for the various Hell features. 
However, the velocities quoted in Table 3 are certainly 
entached of considerable uncertainties (e.g., because 
of blending of TVS subpeaks) making difhcult a clear 
statement as to whether the variability takes place at the 
same characteristic velocities in different lines. 
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1.5 




1002.956 
1342.848 



4580 4600 4620 4640 
Wavelength (A) 



Fig. 5. Superposition for the spectral range encompassing 
Nv AA4604, 4620 of the nightly means of 1995 October 
12 {solid line) and 1996 September 16 {short-dashed line). 
The mean heliocentric Julian dates of the observations (- 
2,449,000) are indicated. 



3.4- The P Cygni Profile Variability 

The TVS analysis also reveals substantial variability at 
the location of the P Cygni absorption component of N v 
A4604 (FigJJ). This mainly results from the transition of 
the N V A4604 feature from a pure emission line-profile in 
1995 October to a P Cygni line-profile in 1996 September 
(Fig.^). No clear variations in the strength of this absorp- 
tion trough on a daily timescale are observed. 

An enhanced peak in the TVS can also be found at 
the location where the carbon triplet C IV A5806 and 



He I A5876 merge (Fig.Q). The large daily changes af- 
fecting Civ A5806 and/or Hei A5876 in October 1995 
and September 1996 are illustrated in Figure ^ (this phe- 
nomenon is not observed in 1996 November). The high 
level of variability observed at this particular location is 
likely due to the superposition of two distinct type of vari- 
ability: (a) red-wing variability of C iv A5806 as observed 
in other spectral features (e.g.. Hen A4686; Fig.^; (b) 
variations in the strength of the P Cygni absorption com- 
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Fig. 6. Superposition for the spectral range encompassing C iv A5806 and He i A5876 of nightly means observed in 
1995 October {left panel), 1996 September {middle panel), and 1996 November {right panel). The mean heliocentric 
Julian dates of the observations (- 2,449,000) are indicated in each panel. 



ponent of He i A5876.0 The latter interpretation is sup- 
ported by the fact that the projected velocity of the peak 
in the TVS {v ^ - 1980 km s~\ referred to the Hei rest 
laboratory wavelength) matches fairly well the wind ter- 
minal velocity (wqo ~ 2135 km s^^). Considering the rel- 
atively modest level of variability affecting the red-wing 
part of the line profiles in WR 1 (Figji]) , it is likely that 
changes in the strength of the P Cygni absorption compo- 
nent of the He i feature mostly contribute to the observed 
changes. 



3.5. Correlated Line-profile Variations in Different 
Lines? 

A possible correlated pattern of variability in two differ- 
ent spectral features has been investigated by calculating 
the Spearman rank- order correlation matrices (see, e.g., 
Johns fc Basri 1995 ), whose elements r{i,j) yield the de- 
gree of correlation between the line intensity variations 
at pixels i and j in each line profile, respectively. In the 
case of perfectly, positively correlated variations, a ma- 
trix unity is obtained. The correlation matrices of He ii 
A5412 with Heii A4686, Hen A4859, and Nv A4945 are 
shown in Figure Q in the form of contour plots, where the 
lowest contour indicates a significant positive (or nega- 
tive) correlation at the 99.5 % confidence level (the other 
spectral lines were not covered enough to be included in 
the analysis). These matrices are displayed in the pro- 
jected velocity frame (referred to the line laboratory rest 
wavelength) . A significant positive correlation is generally 
found between the pattern of variability of Hen A5412 
and He ii A4686. The same conclusion, however restricted 
to the velocity range (- 1000, -I- 1000) km s~^, holds for 



^ Such P Cygni absorption troughs of He I A5876 are com- 
monly (and more clearly) observed in other WN stars (e.g.. 



Robert 1995 



the variations affecting Hen A5412 and Heii A4859. In 
contrast, the (weak) variations of Nv A4945 are appar- 
ently not linked to those of Hen A5412 (the same is true 
for C IV A5806 and He i A5876) . Various factors are suscep- 
tible to mask a potentially correlated pattern (e.g., noise, 
blends). Also, due to the stratified nature of WR winds, it 
is conceivable that (a) two given features present similar, 
yet time-delayed patterns of variability, and (b) the vari- 
ability is mainly restricted to the line formation regions 
of the Hen ions and does not extend where the bulk of 
the N V A4945 emission originates. This would also result 
in an apparent lack of correlation. Overall, the results of 
this analysis tend to show that the Hen features vary in 
a fairly similar fashion. The small number of spectral fea- 
tures included in the analysis prevents, however, to draw 
at this stage more general conclusions. 

3.6. Centroid, FWHM, EW, and Skewness Variations 

We have investigated the time-dependent changes in the 
global line-profile properties by calculating the centroid, 
FWHM, EW, and skewness of Hen A4686, the strongest 
line. The centroid and skewness were calculated as the first 
moment, and the ratio of the third and the (3/2 power of 
the) second central moments of the portion of the profile 
above two in units of the continuum (in order to avoid 
the contribution of blends). The FWHM was determined 
by a Gaussian fit to the entire line profile. The EWs were 
calculated by integrating the line flux in the interval 4643- 
4780 A. One can note that the mean EW of He n A4686 
increased from 313 to 343 A from 1995 October to 1996 
September. This 10 % increase may be explained by an in- 
trinsic EW variability and/or by long-term changes in the 
stellar continuum flux. As can be seen in Figure H {right 
panel), the data obtained in 1996 September generally 
show coherent time-dependent variations. This is espe- 
cially clear for the skewness time series, with the timescale 
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of the variations being of the same order of magnitude as 
the variations of the continuum fiux shown in the upper- 
most part of Figure ^ namely, in the range 5-7 days (note 
that, unlike the EW, the skewness is insensitive to changes 
in the continuum flux level). On the other hand, however, 
and although this may be induced by the paucity of the 
data, no clear time-dependent behavior is noticeable in 
the 1995 October dataset {left panel). 

4. Comparison with Previous Studies 

Niedzielski (1996a) reported on a spectroscopic pattern 
of variability of WR 1 very much reminiscent of the 
one presented by the apparently single WN 5 star WR 
6 that displays phase-locked (although strongly epoch- 
dependent) spectral variations according toV ^ 3.77 days 
(Morel et al. 1998, and references therein). The results of 
our investigation of the spectral variability of WR 1 are 
broadly consistent with this suggestion. In particular, as 
general features of the spectroscopic pattern of variability, 
one can note in both objects the substantial variations of 
the absorption component of the optical P Cygni profiles 
(readily observable in N v AA4604, 4620 where the absorp- 
tion trough occasionally disappears; Fig.||) or the coherent 
time-dependent changes in the global line-profile proper- 
ties (e.g., skewness; Fig.^). Similarities can also be found 
with WR 134, another very rare ex ample of single- line WR 
star displaying cyclical variations ( Morel et al. 1999|) . 

Evidence was also presented by Niedzielski (1996a) for 
cyclical (according to P « 2.667 days), correlated changes 
in the EWs of Hen A4686 and Hen A5412. The data pre- 
sented in the present paper do not, however, support a 
claim of such periodicity. In particular, a period search 
in our EW data (but also in our centroid, FWHM, and 
skewness data; Fig.B|) yields no significant signal at the 



expected frequency. Also, this 2.667 day period can hardly 
account for our global light-curve morphology (Fig.|l]) . In a 
speculative vein, this apparent disagreement between our 
results and those presented by Niedzielski (1996a) might 
point to a multiperiodic nature of the variability in WR 1, 
with the 2.667 day period being present in 1994-1995 but 
being dominated by other longer (non)cyclical processes 
during our observations. It has to be noted that based 
on a more comprehensive analysis of the data, Niedzielski 
(private communication) recently questioned the strictly 
periodic nature of the spectral variations. 

As already noted by Niedzielski (1995, 1996b) and 



Wessolowski & Niedzielski (1996), our data provide evi 



dence for a more gradual spectral pattern of variability 
than for WR 6; a fact which may imply that the period 
(if any) is longer than 3.77 days. In support of this view, 
the centroid, EW, and skewness data of He n A4686 show 
coherent time-dependent patterns of variability (Fig.^, 
with timescales in the range 5-6 days after a formal pe- 
riodicity search. However, we caution the reader that this 
period should not be taken too literally because of the lim- 
ited time sampling of these data. In particular, a cyclical 
pattern in this range, as proposed by Lamontagne (1983) 
and Moffat fc Shara (1986)| , is clearly inconsistent with 
the lack of continuum flux variations observed after HJD 
2,450,353 (Fig.0) unless one invokes a sudden period of 
"quiescence" after this date, with the small amplitude of 
the variations masking any periodic patterns. 

5. On the Presence of a Companion 

Although a cyclical pattern of variability has yet to be un- 
ambiguously established in WR 1, general considerations 
are made below regarding the model assuming an orbiting 
(unseen) companion as the origin of the variability. 
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5.1. A Non-degenerate Companion? 

Some constraints can be set on the mass, M^., of this pu- 
tative companion on the basis of the centroid measure- 
ments of He II A4686 presented in Figure H. Assuming 
that these variations are entirely attributable to orbital 
motion (i.e., Kwr ~ 70 km s~^), one can explore the al- 
lowed values of (V, Mi,). The solid lines in Figure ^ show 
the result of this investigation for three illustrative val- 

o 

ues of the orbital inclination (z = 30 , 60 , and 90 ). In 
these calculations, we a ssume a circular orb it and a mass 
for WR 1 of 9.1 Mq ( iHamann et al. 1995]) . For a wide 

o 

range in orbital inclination (z <; 30 ) and in orbital pe- 
riod {V ^ 20 days), an upper limit for the companion's 
mass of 15 Mq is derived. Assuming the companion to be 
a main sequence star, this constrains the spectral type to 
be lat er than Bl. In this case, the companion's wind is too 
weak ( Grigsby fc Morrison 1995 ) to induce wind- wind col- 
lision effects that may induce the large spectral changes 



observed in WR 1. For a system observed nearly face- 

o 

on (i ^ 30 ), larger masses are evidently consistent with 
the KwR. value adopted above. However, no direct evi- 
dences (i.e., photospheric lines in the integrated spectrum 
or dilution of the WR continuum) support the presence 
of a luminous, early-type companion. In the presence of a 
non-degenerate companion, one would also expect a flat 
or eclipsing-type light curve. Yet, the opposite behavior 
is observed in Figure |^. These arguments strongly argue 
against the presence of a non-degenerate star orbiting WR 
1. 



5.2. A Collapsed Companion? 

For a canonical mass of the companion as a neutron star. 
Mi, « 1.4 Mq, improbably small periods below one day 
are consistent with the adopted Kwr value (Fig.^) . How- 
ever, it has to be kept in mind that the centroid measure- 
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ments of He ii A4686 included the highly variable upper- 
most part of the profile which is unlikely to purely reflect 
orbital motion. Therefore, this value is probably 

grossly overestimated. If one considers lower values 
(as Lamontagne 1983: Kwr = 22 ± 5 km s~^; the allowed 
values of [V, M^] are shown in this case by short-dashed 
lines in FigJ^), much larger periods (as in fact suggested 
by our study) are allowed. Thus, these considerations are 
not sufficient by themselves to rule out the presence of 
a neutron star companion (a black hole companion, on 
the other hand, would require long periods and/or small 
orbital inclinations). ^ 

In the presence of a collapsed companion, one may 
expect fairly strong, accretion-type X-ray emission. Ear- 
liest observations by the HE AO A-1 experiment gave an 
upper limit Lx ~ 4.4 x 10^3 g^g g- i on the emission 



in the 0.5-20 keV range (|Helfand 19801) . A value Lx = 
7.07 ± 2.85 X 10^2 gi-g g-i in the 0.2-2.4 keV range has 
been reported for WR 1 during the ROSAT all-sky survey 



* Considering a lower Kwr value in §5.1 does not qualita- 
tively modify our conclusions. 



( [Pollock et al. 1995| ). WR 1 is a fairly strong X-ray emitter 
compared to other (apparently) single WN stars. However, 
its emission is by no means unusual when only consider- 



ing the WNE-s subclass (Wessolowski 1996). Two subse- 
quent pointed ROSAT PSPC observations showed that 
a satisfactory fit to the X-ray spectrum can be achieved, 
either with a Raymond- Smith thermal plasma of about 
1 keV or by the model developed by Baum et al. (1992)| 
assuming a mixture of "cool" (i n radiative equilibrium) 



and "hot" (shock ed) material ( Wessolowski et al. 1995 



Wessolowski 1996 ). This picture is consistent with our cur- 
rent understanding of the X-ray production in bona fide 
single WR st ars, as being due to rad iatively-induced insta- 
bihties (e.g., [Willis fc Stevens 1996]) . This "normal" level 
of X-ray emission from WR 1 does not constitute, how- 
ever, a decisive argument against the presence of a col- 
lapsed companion, as the accretion of the wind material 
onto the neutron star is known to be inhibited in some 



X-ray binaries (e.g., Zhang et al. 1998). 

Because of the spiral-in process that massive close bi- 
naries are believed to experience in the course of their 
evolution, one is led to expect periods of some hours for 
systems made up of a WR star and a compact companion, 
not days (e.g., De Donder et al. 1997). A period of 4.8 hr 
is observed in Cygnus X-3, the prime candidate for a WR 



compact companion system (van Kerkwijk et al. 1996) 



6. Concluding Remarks 

The single-star hypothesis appears appealing when con- 
sidering the similarities between the optical spectral pat- 
tern of variability of WR 1 and the ones of the two 
peculiar stars WR 6 and WR 134. In this respect, al- 
though the light-curve morphology of WR 1 (a "bump" 
followed by a plateau) has, to our knowledge, no ex- 
ample among the WR populati o n (Moffat & Shara 



1986; [Lamontagne fc Moffat 1987[ ; [Antokhin et al. 1995 
Marchenko ct al. 1998b|) , such a well-defined light-curve 



pattern is reminiscent of what is observed in WR 6 
( Robert et al. 199^ Such a repeatable pattern would, 
however, imply an unlikely large value for the period in 
the context of this single-star hypothesis (V <; 18 days). 

Since the substantial depolarization of the emis- 
sion lines observed in WR 6 and WR 134 is gener- 
ally taken as evidence for an equatorially-enhanced out- 
flow (Schulte-Ladbeck et al. 1991, 1992), revealing the 



same phenomenon in WR 1 may lead to the interesting 
suggestion that the occurence of large-scale line-profile 
and photometric variations (and thus possibly of az- 
imuthally structured outflows) in single WR stars is some- 
how linked to the existence of a wind-compressed zone 



(Ignace et al. 1996) 



^ Curiously, a very similar light-curve morphology (both 
in terms of the timescales involved and of the amplitude of 
the variations) has been noticed in the Be star FV CMa 



(Balona et al. 1992) 
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